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A spin- 1/2 Anderson impurity in a semiconductor quantum well with Rashba and Dresselhaus 
spin-orbit couplings is studied by using a variational wave function method. The local magnetic 
moment is found to be quenched at low temperatures. The spin-spin correlations of the impurity and 
the conduction electron density show anisotropy in both spatial and spin spaces, which interpolates 
the Kondo spin screenings of a conventional metal and of a surface of three-dimensional topological 
insulators. 



A spin-1/2 impurity with large on-site repulsion and 
occupation energy well below the Fermi level of a con- 
ducting band will form a local magnetic moment At 
low temperatures, the impurity state hybridizes or cou- 
ples with the conduction electron states and the local 
magnetic moment is completely screened by a conduc- 
tion electron spin cloud surrounding the impurity, which 
is known as the Kondo effect [2]. To intuitively under- 
stand the Kondo effect, we may consider N-1 conduction 
electron and one impurity electron. In the absence of the 
hybridization, the impurity spin is free to form a local 
moment due to the on-site Coulomb repulsion at the im- 
purity site. In the presence of a small hybridization, the 
virtual process of the hopping between the conduction 
and impurity states favors a spin singlet state of the total 
N electron system. More detailed Wilson's renormaliza- 
tion group analysis shows that the energy gain is of order 
of the Kondo temperature, which has an exponential de- 
pendence to the coupling jSj . After decades of intensive 
study, the Kondo effect as a quantum impurity problem is 
well understood^. The Kondo spin screening cloud has 
long been theoretically predicted. For more recent stud- 
ies, see for exmaple However, direct experimental 
observation of the Kondo screening cloud turns out to 
be more difficult. It was proposed that nuclear magnetic 
resonance could be use to measure the spin polarization 
of conduction electron spins of proper metallic sample in 
dilute magnetic-impurities (Fe for example) by applying 
a small magnetic field. The effect was later found too 
tiny to be detected 0, due to the smaU energy scale or 
the large Kondo spin cloud involved in the screening!^ 
[T]| . Recent development on the spin-resolved scanning 
tunneling microscopy may open a new possibility to de- 
tect the Kondo screening spin cloud in experiments by 
directly measuring the correlation between the impurity 
spin and the conduction electron spin density. This re- 
vives our interest in study of Kondo spin screening cloud. 

In a metal or a semi-conductor, the spin-orbit coupling 
breaks spin SU(2) symmetry. The spin of an electron feels 
an effective momentum-dependent magnetic field and is 
no longer conserved. A natural question is effect of the 
spin orbit coupling to the Kond Problem. Recently, the 



Kondo screening of an Anderson impurity in a helical 
metal has been studied [H, The problem is of inter- 
est in connection to the 3-dimensional (3D) topological 
insulator, whose surface states are described by a helical 
metal with Dirac dispersion. The magnetic doping on 
the surface of the 3D topological insulator Bi2Sez has 
recently been realized. Similar to the Kondo prob- 
lem in a conventional metal, the magnetic moment in 
a helical metal is found to be completely quenched at 
low temperatures. However, the texture of the spin cor- 
relations between the conduction electron and the im- 
purity become more complex. Zitko [Tsj was able to 
map the Anderson impurity model in a helical metal to 
a conventional Kondo problem, in which the conduction 
electrons carry pseudo-spins. His work, however, does 
not provide direct information about the spin screening, 
since the pseudo-spin is a complicated combination of 
the real spins. The helical metal may be considered as a 
strong spin-orbit coupling limit of a 2D electron gas with 
Rashba spin-orbit coupling. In a Rashba model, the ra- 
tio of the spin-orbit coupling to the kinetic energy of the 
conduction electron is a tunable parameter. Therefore, 
the model allows us to study the Kondo spin screening 
continuously from a conventional to a helical metal. 

In this paper, we consider an Anderson impurity in 
a 2D electron gas with Rashba and Dresselhaus spin- 
orbit couplings, which are common in narrow-gap semi- 
conductor quantum wells [l^ [l^. There have been a lot 
of activities recently on the Rashba systems in connec- 
tion with the semi-conductor spintronics. The study of 
the Kondo screening in these systems may be of inter- 
est in manipulation of the spins. We use a variational 
wave function to study the Kondo screening. Similar to 
the Anderson impurity models in conventional metal or 
in a helical metal, we find that the magnetic moment is 
also completely quenched in a semiconductor of Rashba 
coupling at low temperatures. The spin-spin correlations 
between the impurity and the conduction electron den- 
sity are anisotropic in both spatial and spin spaces. In 
addition to the diagonal components of the spin-spin cor- 
relations, the off-diagonal spin-spin correlations have fi- 
nite values due to the broken SU(2) symmetry. 
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The system we study is described by a 2D Hamiltonian 
in a; — y plane. In Nambo's spinor representation, it reads 
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cl[a{cr^ky - (Tyk^) + I3{axkx - ayky)]c\^ 
Hh = ^V]^cld + h.c. 

k 

Hd = (ea- ^i)dU+Ud\d^dld^. 

In the above expressions, k = {kx,ky), and He is the 
Hamiltonian for the conduction electrons of effective 
mass m* . a and /? describe the strengths of the Rashba 
and Dresselhaus type spin-orbit couplings, respectively. 
fj, is the chemical potential, cj^. — (cj^.-^, Cjjj.^) is the cre- 
ation operator in spinor representation, and ay are 
the Pauli matrices. H^ is the Hamiltonian for the impu- 
rity states. The impurity is placed at the origin of the 
x-y plane, and the energy level of the impurity state is 
ed, and d^ — dj). The impurity state hybridizes with 
the conduction electrons through H^- 

Let us first focus on the conduction electrons. The 
strength of the Dresselhaus spin-orbit coupling is deter- 
mined by the asymmetric atomic field of the crystal lat- 
tice and it is absent in the crystal with structural inver- 
sion symmetry. The strength of the Rashba spin-orbit 
coupling in a quantum well structure is allowed to be 
tuned by an external gate voltage and it is usually one 
order in magnitude smaller than the kinetic energy of the 
electron gas. 

After diagonalizing He, we obtain the single electron 
energy eigenvalue eks and the corresponding quasiparticle 
operators 7k, 
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where s = ± is the index of the energy band, k ~ 
(k\ + ky)^ and 0k is the angle of the vector {ak^ + 
Pky,aky + Pkx) with the kx— axis. The ground state 
of many electron system H^ in the absence of Hd is then 
simply given by, |*o) = IlfkijGn 7k±|0)' w'^ere |0) is the 
vacuum, and Q denotes the Fermi sea. 

Now we consider the system with the impurity. For 
an isolated impurity, if < /i < + ?7 , the impu- 
rity is singly occupied and a local magnetic moment is 
formed. Due to the presence of Hh, the conduction elec- 
tron hybridizes the impurity state and tends to screen 
the impurity moment at low temperature. To study 
the ground state properties of the system, we adopt the 
trial wave- function method, which was first introduced by 
Varma and Yafet in 1976[13| • The trial wavefucntion cap- 
tures many of the important aspects in Kondo problem, 
and has been widely used in literature. We believe the 



method should at least qualitative or semi-quantitatively 
accurate to study the Kondo spin screening cloud for our 
purpose. The trail wavefunction is given |12|. 
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where k = {k, s} and d„ = -^{e''^k/2^^ _|. ise^'^^/^d;). 
The large U limit is taken so that the double occupa- 
tion of the impurity electrons is completely excluded. 
The parameter oq and are to be determined using 
the variational method which optimizes the ground state 
energy E = {9\H\-^) / {^\^) . Note that the variational 
parameter n depends on both momentum and spin. The 
spin-dependence is via the band index s contained in k, 
which is a function of the momentum. The construction 
of the trial wave- function ^ is similar to that in Ref. 
for a helical metal. In the limit a = /3 = 0, the trial 
wavefunction is reduced to that for Anderson impurity in 
conventional metal. The wavefunction contains the low- 
est energy sector of the single-particle components and 
is expected to grasp the essential physics of the Kondo 
problem. If all kinds of the particle-hole excitations of 
the conduction electrons are included, the wave-function 
would be the exact eigenstate of the system. 

The variational procedure leads to the expression for 
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This will be used in the calculation of the spin correla- 
tions between impurity and conduction electrons. The 
equation to determine the binding energy, defined by 
the energy difference between the state without the hy- 
bridization term Hh and the variational state, namely 



en ' 



ed- fi- E is given by, 
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Afc > means that the variational state has lower energy 
and is stable against the decoupled state. 

In the case of /x » 2p_ (q,2 + p^)^ we can analytically 
solve Equation ([S]) and find 
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where we have assumed a momentum-independent hy- 
bridization strength Vu_ — V, namely we have assumed 
the hybridization is on-site in the real space, a finite bind- 
ing energy shows the independent impurity plus conduc- 
tion electron system is unstable against the hybridization 
of the impurity and conduction electron with opposite 
spins described in Eq. (3). Since the ground state is 
singlet, the local magnetic moment is quenched. 

To investigate the role of the spin-orbit couplings to 
the Kondo effect, we examine the magnetic properties of 
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the system. The central physical quantity we shall study 
is the correlation function of the impurity spin Sd(0) — 
^(ficrd and the spin density of the conduction electrons 
Sc(r) = ict(r)<Tc(r), 



J™(r)^-^(vI/|5«(r)5,"(0)|M/) 
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that there is also the orbital part of the Kondo screening, 
which contributes to the screening of the impurity spin. 
In the limit the Rashba coupling (assumed to be positive 
without loss of generality) is small, we have analytic ex- 
pressions below for these correlations (hereafter we use 
the units given in the caption of Fig. [T] for brevity) , 



where Af — X^Kgsi "^k possibility for there is a spin 

on the impurity and u,v = x,y,z are the spin indices. 
Juv have the same expression as in the Ref. We 
write them explicitly down here, 

J..(r) - -i|S(r)p + i|^(r)p 
J..(r) = A\Bir)\^-^ReA\r) 
Jyyir) = -i|S(r)p + ii?e^2(^) 

Jxy{r) = Jyx = -^ImA'^ir) 
Jxzir) - -J,x - i/TO(^(r)S(r)) 
Jyzir) = -J,, = -ii?e(^(r)S(r)) (8) 
In the above equations, 

^(r) = ^se*(''-^+'''')a«, 



Firstly, we study the system with a pure Rashba cou- 
pling {a 7^ 0, (3 = 0). From the expression for A and 
equation it is easy to check that ^(r) = in the 
limit a = 0. In Figll] we plot the spatial distribution 
of all the spin correlations between the impurity and 
the conduction electrons. We see that the Jzz is always 
isotropic about the origin point while its density around 
this point decreases as a increases. For a ^ 0, Jxx and 
Jyy are anisotropic and the anisotropy becomes stronger 
at a larger a, as we can see from Fig. 2. In sharp con- 
trast with the system without the spin-orbit coupling, 
all of the off-diagonal spin correlations Juv where u ^ v 
are nonzero and they become larger for larger a. These 
results are closely related to the spin SU(2) symmetry 
broken. 

Since the total spin of the conduction electrons is not 
a good quantum number, the description of the spin sin- 
glet of the total spin of the conduction electrons and the 
impurity spin is no longer appropriate here. A proper 
quantity to measure the spin part of the Kondo screen- 
ing in the present systems is then the correlation be- 
tween the impurity spin and the total spin of the con- 
duction electrons, namely Iq = X)u ^u, where u = x,y, z 
and lu = J dvJuu{Y)- ReA^{Y) has a d-wave symmetry 
in space and it has no contribution to Ix and ly. We 
have Ix = ly and usually they are not equal to Iz- Note 
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From the above equations, at a = 0, we have Ix = ly = 
Iz — — J and Iq = — |, which is the limiting case for the 
Kondo screening in conventional metal. For more general 
values of the spin-orbit coupling, we show our numerical 
results of the correlations in Fig. 2. As we can see, Iz, Ix 
or ly increases from —1/4, initially linearly at small a, 
then to a saturate value of = or /j; = = — 1/8 at 
large a. Note that the value of a may be tunable in ex- 
periments by applying a gate voltage. We remark that a 
small increase of a may suppress the spin correlation sig- 
nificantly and this behavior is enhanced with smaller A;, 
as indicated in Eqn. ©. Note that Iz = implies that 
z-component of the impurity spin is unscreened globally 
by the spins of the conduction electrons (As discussed 
in [T2j , the z-component of the impurity spin is actually 
completely screened by the z-component of the orbital 
angular momentum of the conduction electron). The 
large a limit of /q is — |; which is exactly the value for the 
Anderson impurity in the helical metal with the chem- 
ical potential below the Dirac point. Here we connect 
the Rashba system with the helical metal. The helical 
metal exists on the surface of a 3-dimensional topologi- 
cal insulator [l9l[20l|. It is the large a (or large m*) limit 
of the Rashba system. For the helical metal, we have to 
take an ultraviolet cutoff which is usually taken as the 
half bulk gap. While for the Rashba system, to choose 
an energy cutoff is not severe, the dispersion is bound 
below and the width of the conduction band is a natural 
energy cutoff. 

Now we study the spin-spin correlation functions in 
the presence of a pure Dresselhaus coupling. The prob- 
lem can be mapped onto pure Rashba coupling case by 
a unitary transformation as examined by Shen et al. 12 1| . 

Under the unitary transformation — ^{ux — (Ty), the 
Pauli matrixes are transformed as ax ^ —Cy, ay ^ —ax 
and az — > —az- Therefore under such a transformation 
the Rashba and Dresselhaus couplings in the Hamilto- 
nian ([T]) are interchanged [21]. If the system possesses 
pure Dresselhaus coupling, it can be mapped onto the 
system with pure Rashba coupling under the action of W. 
Denote the spin correlations of the Dresselhaus systems 
as Juv, then they are related to the ones of the Rashba 
system as j;^(r) = J^^(r), ^^(r) = Jyy{r),J'yy{r) = 
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Figure 1: (color online) Spatial distribution of the spin correlations between impurity and conduction electron density for 
systems with only Rashba spin-orbit coupling (/3 = 0) (al)-(fl) are for Jzz, J^^a;, Jyy, J^y, Jxz and Jy^ with a = 0.3; (a2)-(f2) 
are for the corresponding correlations with a = 0.6. The parameters are A;, = 0.02 and fj, — 1. The length unit is kg^ = 10~^ 

meter, the energy unit is = 3.6 x lO'^eV and a is in the unit of 4^.p^ 
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Figure 2: (color online) Correlation functions, given in Eqn. 
(9) , between impurity spin and total spin of conduction elec- 
trons of the ground state in Hamiltonian (1) with pure Rashba 
coupling. The solid (dashed) lines are for Jz, Ix{Iy) and Iq 
with At = 0.02 ( At = 0.005), and /i = 1. The energy and 
length units are the same as in Fig. [T] 



Jxx{r),J'^y{r) = J^y{r),J'^^{r) = Jy,_{r) and J^Jr) = 

In the presence of both the Rashba and Dresselhaus 
couplings, the energy dispersion is no longer isotropic in 
momentum space as in the Rashba system and so does 
the a^. However, they still possesses reflection symme- 
tries in momentum space about the lines i ky = 0. 
An immediate result is that Jzz{j^) loses its isotropy and 
is symmetric about x ± y — in real space. As shown 
in Fig. [21 for systems with the Dresselhaus coupling, the 
patterns of Jxx, Jyy, Jxz and Jyz are rotated and Jxy 
is remarkably redistributed. At a = /?, the model has 
an interchange symmetry between x and y. In this case, 
Jxx — Jyy and Jxz — Jyz- 



In summary, we have examined a spin- 1/2 Anderson 
impurity in a two-dimensional electron system with spin- 
orbit couplings. The binding energy and the spin correla- 
tion functions are calculated using a trial wave-function 
method, which is widely used to study Kondo problem. 
Being consistent with former studies |23l. [2^ . the mag- 
netic moment of the impurity is found fully quenched at 
low temperatures, as the same as in the conventional An- 
derson impurity problem. However, because of the spin- 
orbit couplings, the spin correlations between the impu- 
rity and the conduction electrons are different from the 
ones in conventional metal. The diagonal components 
Jxx and Jyy are anisotropic in space. The off-diagonal 
components are nonzero and have particular spatial dis- 
tributions. The introduction of the Dresselhaus coupling 
changes the patterns of the correlation functions dramat- 
ically. 

Our theory is variational in nature. The variational 
wavefunction method has been widely used to study 
Kondo impurity problem and the basic results are con- 
sistent with more accurate methods such as numerical 
renormalization group technique. We note that spatial 
extension of the Kondo cloud in the Anderson impu- 
rity model in conventional metal was recently studied by 
Bergmannp^ by using Friedel artificially inserted res- 
onance method, where the screening length is quanti- 
tatively calculated. Our focus in the present paper is 
on the anisotropic spin correlation functions due to the 
Rashba type spin-orbit coupling. We expect that the spa- 
tial spin screening length in our variational calculation is 
also associated with the Kondo temperature, similar to 
Ref . ^3 . In the parameter region we study, the screening 
decay is rather rapid as we showed in the case of the An- 
derson impurity in a helical metal[l3|. We believe that 
our results presented here on the anisotropic screening 
are qualitatively or semi-quantitative correct. Although 



5 



-0.2700 -0.1260 0.01600 -0.2700 -0.1260 0.01800 -0.2700 -0.1260 0.01G0O -0.02000 0.02000 -0.09400 0.09400 -0.09400 0.09400 



4 



• 


• 


• 








A 


(a1) 


(b1) 


(Cl) 


{d1) 


(el) 


(fl) 


4 

-4 


• 

(a2) 


• 

(b2) 


• 

(c2) 


(d2) 


« 

(62) 


(fZ) 




-4 4 


-4 4 


-4 4 


-4 4 


-4 4 


-4 4 



X X X X X X 



Figure 3: (color online) Spatial distribution of the spin correlations for systems with both Rashba and Dresselhaus spin-orbit 
couplings, (al)-(fl) are for Jzz, Jxx, Jyy, Jxy, Jxz and Jyz with a — 0.3 and /3 = 0.1; (a2)-(f2) are for the corresponding ones 
with a — P = 0.3. The other parameters are At, = 0.02 and ^ = 1. The energy and length units are the same as in Fig. [1] 



more sophisticated methods may be needed to provide 
precisely quantitative results of the Kondo spin cloud, it 
may be left for future work. Our study may be of in- 
terest to semi-conductor spintronics, where the Rashba 
system has played important roles. The impurity spin 
may be used to control or to manipulate the conduction 
electrons and the physics of Kondo screening in these 
systems should be interesting. Direct experimental ob- 
servation remains a major challenge. Nuclear magnetic 
resonance does not seem to be promising for the tiny sig- 
nal in large Kondo spin cloud. Recently developed spin 
resolved STM may offer a new route. By adding a small 
magnetic field around the impurity site, the STM may 



in principle measure the spin-dependent local density of 
states to probe both diagonal and off-diagonal spin corre- 
lations. However, it may be still challenging to probe the 
electrons in the quantum well structure, which may re- 
quire long wave optical probe. Nevertheless, the modern 
technique is developing very rapidly, and the interesting 
phenomenon of the Kondo spin cloud, and the anisotropic 
spin screening examined in the present paper will be ob- 
served in the near future. 
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